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ABSTRACT 
We present a detailed dynamical analysis of the thin, extended stellar stream that encircles the 


disk galaxy NGC 1097. 


Within a fully statistical framework, we model the surface pgpness map 
of the so called ‘dog leg stream’, using mock streams generated as in pose) (2015 (2015) 
deep imaging data obtained with the CHART32 telescope for the Stellar Tida 


, and new 
tream ee We 


reconstruct the post-infall evolution of the progenitor galaxy, which, we find, has experienced three 
pericentric passages, and lost more than two orders of magnitude in mass. At infall, 5.4 + 0.6 Gyr 
ago, the progenitor was a disky dwarf with a mass of log, )(m/Mo) = 10.35 + 0.25 within its tidal 
radius, r, = 3.441 kpc. We illustrate how the sharp 90 degree turn in the stream, identifying the ‘dog 
leg’, is the signature of the progenitor’s internal rotation, inclined with respect to the orbital plane 
and imprinted in the escape conditions of the stream members. Today, the remnant is observed to 
be a nucleated dwarf galaxy, with a line-of-sight veley. of Up ie = —30+ 30 kms! with respect to 


los —_ 
NGC 1097, and a total luminosity of 3.3 x signees ea (Galfanni et al J2010). Our independent analysis 


predicts a line-of-sight velocity of up jos = 
logig(m/Mo) = 


743". implying that the Bee nucleus (L 6.9 x 10°Ly.o, |Galianni et al. 
2010) is soon destined to become a low-luminosity Ultra Compact Dwarf galaxy. The progenitor’s 


kms~* and measures its current total mass at 


orbit is strongly radial, with a pericenter of a few kpc, and an apocenter reaching rapo = 150+ 15 
kpc, more than a half of the inferred virial radius of the host, ragq9 = 250 + 20 kpc. We find — 
NGC 1097 has a mass of Mao9 = 1.81?) x 10!? Mo, and its inforted concentration c299 = 6.7773 


is in very good agreement with the expectation of the cold dark matter model. 


We can deseube 


the stream almost down to the noise using a spherical host potential, and find that this would not 


be possible in case of a halo that is substantially triaxial at large radii. 


In the morphology of the 


stream, we can see the logarithmic density slope of the total density profile, y, bending from its inner 
value, (peri) = 1.5 £0.15, and steepening towards large radii. For the first time on an individual 


extragalactic halo, we measure the outer slope of the density profile, 


(0.67200) = 3.9+0.5. This 


demonstrates the promise of the newborn field of detailed, statistical modelling of extragalactic tidal 


streams. 


Subject headings: galaxies: dwarf — galaxies: 


structure — galaxies: 


kinematics and dynamics — 


galaxies: interactions — galaxies: individual NGC1097 


1. INTRODUCTION 


Recent deep wide-field imaging surveys of stellar tidal 
streams have mainly focused on nearby Milky Way-like 
spiral galaxies that were suspected to contain diffuse- 
light, very faint features in their outer regions based on 
existing data ai large scale surveys, such as POSSII 


1997) and the SDSS (Miskolezi et al. 
2011), or ae as (Martinez-Delgado et al.{2008} 
The morphology of the detected tidal streams 
include great circle streams resembling the Sagittarius 
stream around our Galaxy, isolated shells, giant clouds 
of debris floating within galactic halos, spike-like fea- 
tures emerging from galactic disks, and large-scale diffuse 
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structures that may be related to the remnants of an- 


cient, already thoroughly disrupted satellites (Martinez- 
Delgado et al.J2010). This extraordinary variety of mor- 


phologies represents one of the first comprehensive pieces 


of evidence for the hierarchical formation scenarios bie. 
dicted b ee a e.g.,|White & Frenk 


2008} [Cooper 


tellar 


| 
ty 


out the first systematic survey of stellar tidal streams to 
a surface brightness sensitivity of ,-=28.5 mag arcsec”? 
using a network of small, robotic telescopes placed on 
different continents. 

These debris contain precious information on both (i) 
the evolutionary struggles of the progenitor dwarf galax- 
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ies, transformed and partially destroyed by the tides; (ii) 
the formation history of the massive host, its halo popu- 
lations and the structure of its dark matter halo. How- 
ever, rewinding the movie of these accretion events has so 
far been difficult. It is often possible to identify a set of 
properties of progenitor and host that provides a qual- 
itative resemblance to the observed stream 


Delgado et al.|/2008} |2009 
al-[2014}[Martinez-Delgado et al.[2014), 


ration of the wide parameter space is exceedingly expen- 
sive to perform with N-body simulations, meaning that 
the extraction of reliable measurements is not possible. 
Therefore, it remains extremely difficult to provide ro- 
bust confidence limits on several interesting quantities, 
such as the progenitor’s initial mass, infall time and or- 
bital properties, as well as host halo mass, concentration 
and detailed density profile. 

Extended and reasonably cold streams represent one of 
the least degenerate probes of the gravitational potential 
in which they orbit, as their global morphology directly 
traces the density slope and geometrical structure of the 
host galaxy. While it is not correct to assume that a 


stream traces an orbit (Sanders & Binney|/2013) as dif- 
are direc 
Johnston et al./2001 


ferent positions along its leng associated to 
sgradionh mor ital encizy 
it is qualitatively correct to say that, as for 
orbits, different density profiles imply substantially dif- 
ferent shapes. For example, the opening angle between 
successive apocenters is certainly 
a function of the underlying gravitational potential. Fur- 
thermore, in addition to what is easily understood within 
the single orbit interpretation, streams also provide addi- 
tional constraining power, in the shifting galactocentric 
distances of successive apocenters (increasing along the 
trailing tail and shrinking along the leading tail), as well 
as the links between the stream’s total length and width 
and the inner density slope of the host (Anoriscof201). 

The interplay between the central dark matter distri- 
bution and the lifecycle of baryons is one of the crucial in- 
gredients of the galaxy formation recipe. The predictions 
of dark-matter-only cosmological simulations are quite 
clear in this respect: (see e.g., [Navarro e al 10072009), 
but we still lack a thorough understanding of the equi- 
librium between competing processes such as adiabatic 


contraction (Blumenthal et_al./1986;|Gnedin et_al./2004) 
and stellar feedback (Navarro et al.]1996} 
[al.|[2006 Pontzen & Governato|2012) across the galactic 


mass scale. erefore, the inner total and dark matter 
density slopes are currently subject to close observational 
attention, mainly throug 


matics (eg. Cappellari et al [2013] 
erences therein) and strong lensing analyses (e.g 
, 
onversely, it is considerably more dif 
strain the outer density profile of galactic dark matter 
haloes. Deep HI rotation curves of disk galaxies extend 
to ¥ 10 disk scalengths, or ~ 50 kpc (e.g. 
(2008), just beginning to explore the regime r = ro, where 
ro is the scale radius of the halo. There’s also good hope 
for dynamical analyses based on discrete tracers in the 


haloes of massive galaxies: using Planetary Nebulae and 
Globular Clusters it is possible to probe scales of = 100 


kpc, henceforth reaching a fraction of the virial radius 
(Napolitano et al.|/2011} |Pota et al.|/2013} |Brodie et al. 


Fic. 1.— Luminance filter image of the “dog leg stream” around 
NGC 1097 obtained with the CHART32 amateur telescope (see 
Sec. 2.1). The image covers a total field of view of 25.2 arcmin 
x 31.2 arcmin, that at the distance of 16.96 Mpc corresponds a 
projected area of 127 x 157 kpc. The thin white cross identifies 


the remnant of the progenitor of the stream, currently a nucleated 
dwarf galaxy (Galianni et al./2010). 


2014} |Longobardi et al.|/2015). However, modelling has 
to struggle against the intrinsic degeneracies of pressure 
supported equilibria (e.g. 1982), and 


while resorting to the joint modelling of distinct subpop- 


ulations is helpful in inferring the properties of the inner 
profile (Agnello eal 2014) Napolitano et a 2014, mea 
suring the mass and density profile at large radii might 
just be too difficult. 

As a consequence, observational constraints on the 
outer density profile of galaxy-sized dark matter haloes 
are extremely scarce. (2007) have used a 
combined strong and weak lensing analysis on a sample 
of 22 early type massive lenses (aperture-averaged cen- 
tral velocity dispersions between 200 and 300 kms~*), 
and inferred the average dark matter distribution out to 
= 300 kpc, or © 0.75ryi;. They find that the weak lensing 
signal can be satisfactorily described by a density profile 
having an outer slope of y, = —3. However, the radial 
range actually probed does not extend much beyond the 
regime ypm(r) * —2, nor the available data is sufficient 
to allow for a comparison between models with different 
asymptotic density slopes. Analogous weak lensing anal- 
yses on stacked data are also available for more massive 
systems, as galaxy groups and clusters. 
(2006) sl Ut ta 2011 can see the total den- 
sity slope bending from an inner value compatible with a 
typical NFW profile, y; = —1, and steepening in a man- 
ner that can be well described by a profile with y, = —3. 
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Differently from the constraints we just mentioned, 
which use stacked data and only reveal the average prop- 
erties of haloes, extended tidal features offer us the op- 
portunity to explore this dark matter dominated realm 
in individual galaxies, without having to worry about 
the additional biases introduced by stacking and by the 
homology hypothesis. Detailed modelling of extragalac- 
tic tidal features is definitely a field in its infancy, but 
the quick developments in the description and modelling 
of tidal streams (see e.g., 
and references therein), together with the observa- 
tional race towards such low surface brightness targets 
2014} |Duc et al.//2015) ensure fast and promising devel- 
opments. Expectations also include a clearer picture of 
the evolution of dwarf galaxies, by the reconstruction of 
the post infall evolution of the progenitors of the studied 
streams. 

The extended jet-like features apparently emanating 
from the central regions of NGC 1097 are not a recent 
discovery (e.g. Asp] 1970) however, their real nature has 
remained elusive for almost two decades, until their pos- 
sible associations with radio jets was finally ruled out 
(Wehrle et.al. and references therein). The four 
plumes in the X-shaped stream (see Fig.|1) reach out to 
projected distances of over 60 kpc, extending well beyond 
the barred disk of NGC 1097. The apparently longer 
plume displays a peculiar, sharp 90 degree turn, which 
has originated the name of “dog leg stream” 
fal. [2010). The same authors identify a nucleated dwar 
— total luminosity ~ 3.3 x 10’Ly.o and half light ra- 
dius of ~ 300 pc — well aligned along the ‘leg’ of the 
stream, as indicated in Fig. They spectroscopically 
confirm its association with NGC 1097, and show that 
the stream appears to be emanating directly from this 
“dwarf Spheroidal-like” object, which is then recognised 
as the remnant of the progenitor. 


The minor merger scenario has been explored by 
don & Wallin| (2003), which try to roughly reproduce 
t 


e stream’s morphology using N-body simulations of the 
disruption of a disky dwarf. They conclude that the pe- 
culiar X-shape of the stream in NGC 1097 is due to the 
strong disk component of the host, as they do not seem 
to find viable models when using spherical potentials. 

In this paper, we use new deep imaging data of the 
stream in NGC 1097 and a new dynamical technique to 
reconstruct this accretion event within a fully statisti- 
cal framework. We model the surface brightness of the 
stream using mock streams orbiting static, spherical po- 
tentials (Amore 2015) These are generated within 
a mixed parametric/non-parametric approach, which al- 
lows us to follow the mass evolution and disruption of 
the progenitor. This technique is computationally effi- 
cient enough to allow for a maximum likelihood imple- 
mentation and a systematic exploration of the param- 
eter space, so that reliable confidence intervals can be 
inferred for the properties of both progenitor and host. 
Section 2 presents the available data; Section 3 provides 
an overview of the generative technique; Section 4 shows 
that the projected shape of an orbit is sufficient to con- 
strain the underlying gravitational potential; Section 5 
and 6 present the results of this paper, respectively per- 
taining host galaxy and progenitor dwarf galaxy; Sec- 


Fic. 2.— The masks used to construct the surface brightness 
map of the stream. Left: the host galaxy, foreground stars and 
bright background galaxies are masked out. Right: the region of 
the stream is isolated to learn the properties of the noise in the 
remainder of the image. 


tion 7 comments on the limits of this modelling approach 
and discusses our results. 


2. OBSERVATIONAL DATA 
2.1. Observations 


Deep imaging of the field around the spiral 
galaxy NGC1097 was collected at the remote obser- 
vatory CHART32?} (Cerro Tololo Interamerican Ob- 
servatory,Chile) with a 80-cm aperture f/7 corrected 
Cassegrain telescope as part of the Stellar Tidal Stream 
Survey (PI. Martinez-Delgado). A FLI PL-16803 CCD 
camera was used and provided a pixel scale of 0.331! /px- 
1 over a 22’ x 22’ field of view. Two overlapping sets of 20 
x 1200 s individual image frames were obtained remotely 
with a Baader luminance filter over several photometric 
nights between September and November 2013. Each in- 
dividual exposure was reduced following standard image 
processing procedures for dark subtraction, bias correc- 
tion, and flat fielding adopted for the larger stream sur- 
oe (Martney-Delgado et a 20T0) ate 
combined to create a final co-added image with a total 
exposure time of 48,000 s. The resulting luminance im- 
age is displayed in Fig. 1. A detailed description of the 
photometric calibration of the data, together with our 
near-infrared and optical photometric follow-up observa- 
tions with 8-meter class telescope of this stream, will be 
presented in a companion paper (Martinez-Delgado et al. 
in preparation). 


2.2. Extracting data and uncertainties 


We use the data in the luminance filter to extract a 
map of the surface brightness of the stream and of the 
associated uncertainties. Note that we do not need a 
calibrated map of the absolute luminosity of the stream, 
but only a map of the projected locus occupied by the 
streaming material, and of the relative brightness of dif- 
ferent regions along such locus. In order to simplify the 
comparison with models, we construct this map by ras- 
terising the original image, i.e. by dividing the area cov- 
ered by the data in larger ‘effective pixels’. We associate 
an average brightness — with uncertainty — to each of 
these. There are two main difficulties to overcome in do- 
ing this: (i) isolating and excluding all compact bright 
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sources present in the data, such as foreground stars and 
background galaxies; (ii) estimating and accounting for 
the smooth variations in the non-flat background lumi- 
nosity. 

We solve issue (i) using a twofold strategy. First we 
take much care in manually masking out all foreground 
stars and background galaxies brighter than some thresh- 
old. The result of this operation is illustrated in the left 
panel of Fig. [2] which shows the resulting mask (excluded 
region in white). This concurrently masks out the host 
galaxy, the progenitor of the stream and the region in the 
upper left corner of the data. The reason for excluding 
the latter portion of data is that this contains the poorly 
determined apocenters of the two upper plumes of the 
stream. Because of scattered light from a bright fore- 
ground star and the intrinsically low surface brightness 
of this section of the stream, the exact apocenters of these 
plumes are highly uncertain. As a consequence, rather 
than imposing incorrect constraints through possibly in- 
accurate data, we prefer to leave complete freedom as to 
the extension of these plumes, and mask this region out. 
Second, in the original image, we systematically exclude 
all single pixels that are brighter than a given threshold. 
This is considerably higher than values reached within 
the stream, so this final step simply takes care of mask- 
ing all spurious sources that had survived the first mask. 
We refer to the final non masked region using the symbol 
M. 

In order to subtract the background, we smooth it 
through a classical Gaussian kernel. For a given size 
of the effective pixels, we calculate average brightnesses 
and uncertainties using the following classical estimates 


= (u)i —(M)G 4 » (1) 
Nj dv? = ((u")i — (HE) + (Hea — (ui) > (2) 


Yj 


where 1; indicates the brightness of the effective pixel 7 
and dv; its uncertainty, while y is the brightness of pixels 
in the original data; (-); is the average over the area of the 
effective pixel i and (-)gv,; indicates the average centred 
on the effective pixel 7 and obtained using a Gaussian 
kernel — with characteristic size a few times larger than 
the size of the effective pixel itself. Note that both av- 
erages only use pixels in the original image that are not 
masked out, i.e. pixels in M, so that N; is the number 
of such pixels used for the effective pixel 7. 

So to learn about the properties of the noise, we sepa- 
rate the available data in two regions, as from the mask 
in the right panel of Fig. |2| isolating the stream in the 
region M,t;. In the complementary region, masked out 
in the right panel of Fig. [2] we impose that 


Vi 
where p indicates the probability distribution of the 
brightness of the effective pixels, normalised by their 
uncertainty. In other words, we rescale the uncertain- 
ties calculated using eqn. (2) so that the distribution in 
eqn. is as similar as possible to a standardised Gaus- 
sian, as it would be if outside Mgt,, where the stream is 
absent, the image is dominated by noise. In practice we 


140, 


68 
14 42352216167 354121 


pdf 


0 5 10 


y;/ndv; 


Fic. 3.— Upper panel: the distribution of the noise in the region 
where the stream is absent, in grey, together with an histogram 
of the significance of the signal in the region of the stream (see 
Fig.|2). Lower panels: the map of the surface brightness in the re- 
gion of the stream (left), is shown together with the corresponding 
uncertainties (right). 


- (@) (2)] (4) 


and find that, for an effective pixel size that is a fraction 
of the stream’s width, values n ~ 3 are needed. The 
grey shaded area in the upper panel of Fig. [3] shows the 
resulting distribution of the ‘signal’ in the region where 
the stream is absent, masked out in the right panel of 
Fig. This well compares to a Gaussian distribution 
with zero mean and unit standard deviation, also shown 
as a full black line. Therefore, we apply the scaling re- 
quired by eqn. to all effective pixels, including those 
in the region containing the stream, M<gtr. 

The final result of this procedure is shown in the lower 
panels of Fig. using the same grayscale for both av- 
erage brightness and associated uncertainties, these il- 
lustrate the map of the effective pixels in the stream’s 
region M,t;._ The stream’s morphology can be clearly 
recognised, with values 1;/nd1;, reaching up to over 14 
for this particular effective pixel size, corresponding to 
a side of 2 kpc in the sky at the distance of NGC 1097 
(D=16.96 Mpe, [Roribalst etal] 2003), The effective 
pixels characterised by significant uncertainties are due 
to the Poisson noise of eqn. je as they are close to re- 
gions masked out of M, N; decreases considerably. The 
complete distribution in »;/ndv; of the region Mgt; is 
shown in the upper panel of Fig. [3] with a (scaled) red 
histogram. 

Finally, we use the original image to measure the cur- 
rent position of the remnant of the progenitor in the 
plane of the sky. With respect to centre of NGC 1097, 
we find 


(Zp.e1;2p,ep) = (51.5, —33, 8) + (0.3,0.3) kpe, (5) 


use 


Model of the ‘dog leg’ stream 5 


where the directions 6; and é2 are respectively the hor- 
izontal and vertical directions in Fig. |1} in the plane of 
the sky. This corresponds to a projected distance of 


Ry =61.6+40.4 kpc . (6) 


3. STREAM GENERATOR 


3.1. The technique and the free parameters 


We refer to (2015), for brevity A15, for a 
complete description of the method used to generate 
mock tidal streams, and report here a brief summary of 
the strategy and a description of any additional feature. 

As in A15, the underlying total gravitational potential 
is generated by a spherically symmetric broken power 
law profile 

—(Yo-¥)/2 


wo POP. o 


with inner and outer slopes y; and yj, and break radius 
ro. Orbits in such potential are described using a pur- 
posely optimised orbit library, that covers the parameter 
space 


VE pas 

VY € [2.1,5 

ge Ot) (8) 
logio(Te/T0) € [-2, 2] 


where j is the classical circularity and the ratio r./ro 
is used as a convenient substitute to orbital energy, 
r-(E/Gporg) being the radius of the circular orbit with 
the dimensionless energy E/Gpor?. As fully described in 
A15, this technique bypasses all ‘on the fly’ integrations 
of the differential equations of the motion and, as a con- 
sequence, is such that evolution of a set of particles to 
any given time is equally expensive. 

Mock streams are obtained by releasing particles along 
the orbit of the progenitor (rcp»/ro,jp), starting from 
infall (operatively, the first apocenter) and until the 
present time, according to a parametrized shedding his- 
tory Psn(t). This is the sum of Gaussian distributions, 
one for each pericentric passage, with free mean shedding 
time, characteristic spread and shed mass. 

The initial conditions at escape (dr, dv) can also be 
varied at will, so to mimic the mass evolution of the pro- 
genitor or the disruption of a system with internal angu- 
lar momentum. Particles are released in the proximity 
of the instantaneous tidal radius of the progenitor 


Gm 3 
—— (a = 55 58) (9) 


the size of which depends on the progenitor mass m and 
on the instantaneous galactocentric distance. Particles 
leave the progenitor from apertures in the surface of the 
effective potential, the size of which is connected to their 
energy, or equivalently to the escape velocity dv. Esti- 
mates for this link have been provided in A15, and we 
report them here for reference: 


ea ry |Ov|? 

Qp — t G ’ 
oe ee (10) 
V2 = Wey! +1 


where @, and @, are the semiaxes of the elliptical es- 
cape apertures, lying in the plane perpendicular to the 
radial direction 7, crossing r;. ¢ is the angular direc- 
tion within the orbital plane, while 2 is perpendicular to 
the orbital plane. As in A15, these quantities are used 
as characteristic scales of Gaussian probability distribu- 
tions, describing the pdfs of the spatial initial conditions 
at escape (see eqns. (60-62) in A15). 

As to the kinematic escape condition, particles are lost 


with the same instantaneous angular velocity of the pro- 
genitor (Kipper et al-2012) Gibbons eta 2014) and an 
additional velocity component dv. Its cylindrical com- 
ponents — referred to the orbital plane — have a Gaussian 
distribution with adjustable mean and dispersion. Using 
the notation of A15, we find that 0, = 0, = 0, = 9g is 
sufficient for the data we have in hand, so that the spread 
in the kick velocities is described by a single function of 
time, o,. On the other hand, we need flexibility in the 
means of the different components. In particular, we use 
the following parametrisation 


(dug) = 6Vrot C08(Pincl) (11) 
(6vz) = SVzot Sin(dinct) cos[y(t) — dorient] 


where 6V,o; phenomenologically corresponds to internal 
rotation in the progenitor and is a function of time; inc 
indicates the inclination between internal and orbital an- 
gular momentum, and @orient allows us to change the ini- 
tial orientation of the disky progenitor, as illustrated in 
Al15. y(t) is the instantaneous orbital location of the 
progenitor. 

It is worth recalling that we only model the release of 
‘stellar particles’ and our mock streams will be used to 
fit the stellar stream alone. Dark matter is also lost to- 
gether with stars: the dark stream is likely much hotter 
and phase-mixes more rapidly as a result of warmer ini- 
tial conditions, as the halo is more extended and has an 
accordingly warmer kinematical profile within the pro- 
genitor. 

Time dependences are treated in a non parametric way. 
Time dependent free parameters are: (i) m(t), which im- 
plies the size of the tidal radius; (ii) o,(t), which is kept 
independent of the progenitor mass, and fixes the spread 
in the kick velocities at escape; (iii) dVrot(t), which ap- 
proximately fixes the time when any disk component of 
the progenitor is disrupted and lost, and the magnitude 
of such signature in the escape velocities. For m and 
go, we simply use a smooth, spline interpolation between 
the values of these functions at each pericentric passage, 
which are treated as free parameters. We find that this 
is enough for our scopes and the current data, and we 
do not need to increase the number of free parameters 
by using a finer mesh in time. For dV,ot, instead, we 
possibly need to resolve finer time scales, and therefore 
use a different approach. We parametrize OV;or(t) as the 
sum of two Gaussians in time, with free means, spreads 
and maximum rotational velocities. This allows us to 
phenomenologically capture at the same time the quick 
disruption of a fragile disk, and/or the durable effects of 
a rotating compact progenitor. 

In conclusion, the free parameters and free functions 
of the model are listed below. 


e The progenitor’s orbit (r-p/ro,jp) and the time 
spent since infall, i.e. the first apocenter, in units 
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[ro] 


[ro] 


Fic. 4.— A typical dataset used to test the constraining power 
of the projected shape of a single orbit. The large panel illustrates 
the projected orbit itself, while the inset shows its corresponding 
rasterised surface brightness map. 


of the progenitor’s orbital time, ting/T;p. 


e The progenitor’s line-of-sight position with respect 
to the centre of NGC 1097 at the current time, 
Lplos — the projected position of the progenitor 
in the sky is fixed by the data within its uncer- 
tainty, see eqn. (5). A possible rotation of the 
orbital plane, @>;p, around the vector identifying 
the position of the progenitor at the current time; 
in practice, this fixes the inclination of the orbital 
plane with respect to the line-of-sight direction. 


e The host total density profile, parametrized by the 
two parameters 7;, Yo. The characteristic density 
po is not directly a free parameter as data are lim- 
ited to relative surface brightness values only, and 
no kinematic information on the stream is avail- 
able. At the same time, the break radius 7 is only 
apparently a free parameter: once the energy of the 
progenitor’s orbit and the three-dimensional posi- 
tion of the progenitor are chosen, rg is univocally 
determined. 


e The shedding history p,,(t), ie. three free parame- 
ters for each Gaussian pdf per pericentric passage, 
minus one, as the integral of the complete pdf is 
normalised to unity. 


e The mass evolution of the progenitor m(t), i.e. the 
mass itself — or equivalently the tidal radius — at 
each pericentric passage. Note that this is treated 
independently from the shed luminous mass at the 
previous point. 


e The velocity dispersion at escape o,(t), ie. its 
value at each pericentric passage. 


e The ordered kick velocities at escape dV; o4(t), ie. 
three free parameters for each of the two used 
Gaussian functions plus two angles for inclination 
and orientation, inc) and Porient- 


yr) 


pericenter 
f 
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Fic. 5.— The inference on the profile y(r) of the total density, 
recovered from fitting the single orbit displayed in Fig./4| The grey 
shaded area represents the 68% confidence region, the 95% confi- 
dence region is not shaded. The red profile illustrates the target 
profile, while the vertical lines display the locations of pericenter, 
break radius and apocenter. 


3.2. Likelihood and model comparison 


We populate each stream with Net: = 1.25 x 10* parti- 
cles, scale it and and rotate it so that the position of the 
progenitor in the plane of the sky at the current time co- 
incides with the one observed in the data. The number 
of particles in each effective pixel n; is calculated, and 
transformed into a quantity that can be directly com- 
pared to the data. This transformation uses an estimate 
of the surface brightness per particle ypp, which is ob- 
tained from the region of the stream M<¢,, highlighted 


in Fig. 2] 
Upp = » v/ S- Ny - (12) 
M 


str Mstr 


For each model .@, identified by its set of free parame- 
ters, this produces the likelihood 
(vi-Vppni)? 
exp [-3 6v24+v2 nj 
£(m)={{ Ss, (13) 


i Qn (dv? + v2,n4) 


based only on the map of the relative surface brightness 
of the stream and in which the Poisson noise of each 
effective pixel is taken into account together with the 
observational uncertainties. 

Given the granularity of the mock streams, samplings 
of the pdf functions that define a given model .@ pro- 
duce slightly different estimates of the likelihood #(.@); 
in other words, #(./) is subject to variations when cal- 
culated repeatedly. This variability disrupts the natural 
exploration of the parameter space, as it causes walkers 
to get stuck in samplings that are fortuitously performing 
better than the average of that model. As a consequence, 
such variations must be suppressed, and we do this by 
tempering the likelihood. For a given effective pixel size 
and number of particles in the stream Nzt,, we use 


logs lV (LMPIT) —(L(M)/T)2| =1, (14) 


where averages are calculated on repeated samplings of 
the same model .@, which, operatively, is taken to be 
the best fitting model. In this way, oscillations of the 
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Fic. 6.— Corner plots for the parameters recovered fitting the single orbit displayed in Fig. Contours indicate the 68% and 95% 


confidence regions, while red lines identify the target values. 


tempered likelihood Y!/7> are limited an order of mag- 
nitude at most, allowing for a smooth evolution of the 
walkers. For Ngty = 12500 and an effective pixel size of 
4 kpc?, we find that T; ~ 12 is needed. 


4. TESTS ON SINGLE ORBITS 


In order to test whether it is possible to constrain the 
density profile of the host using the projected morphol- 
ogy of the stream alone, we perform tests on single orbits. 
A sample of orbits is chosen with properties comparable 
to those inferred for the progenitor’s orbit, which we il- 
lustrate in Section 5. Such orbits are populated with 
~ 10° particles and the resulting surface brightness in 
the plane of the sky is mapped in effective pixels. In 
order to mimic the characteristic width of the plumes 
of dog leg stream, the size of the effective pixels used 
for these tests is taken to be approximately four times 
larger in area than those used to generate the maps of 
Fig. As shown in Fig. |4| this considerably broadens 
the projected profile of the orbits. 

We fit such mock datasets using single orbits from our 
library and find that target parameters are always well 
recovered. Figure [5] illustrates how the measured radial 
profile of the density slope compares with the generating 
one (in red). The grey shaded area corresponds to the 
68% confidence region, contained between the 16% and 
84% quantiles; the 95% confidence region is also shown. 
Vertical lines illustrate the position of orbital pericen- 
ter and apocenter, together with the break radius of the 
underlying density profile. 

Figure iG illustrates in more detail the 68% and 95% 
confidence regions of the parameters that describe the 
main projected morphology of the orbit, and how these 


compare with the target values: as shown by the red 
lines, all values are recovered within one sigma. Most 
parameters are not degenerate and the only noticeable 
degeneracies are those associated with the profile of y(r). 
The inner and outer slopes, 7; and yo, can be slightly 
varied around the best fitting values together with the 
energy of the orbit r./ro, which, as mentioned earlier, is 
a measure of the the break radius. Such degeneracy man- 
ifests itself ina clear pinch point in the recovered profile 
y(r), in Fig. |5| located where the profile first begins to 
steepen, between pericenter and break radius. 


5. RESULTS: THE HOST GALAXY 


We explore the available parameter space using a suite 
of Metropolis-Hastings walkers. We first perform a series 
of tentative explorations employing likelihoods temper- 
ated with high temperatures, using multiples of T;, de- 
fined in eqn. (4p. These walkers quickly swipe through 
very wide areas of the parameter space and approxi- 
mately identify the useful region. Final runs are per- 
formed using T = Ts. 

We repeat the complete analysis on data obtained us- 
ing two different effective pixel sizes — 4 kpc? and 3 
kpc? — so to test the influence of these on any results 
together with the influence of a shifting binning pat- 
tern on the image. Furthermore, we perform the final 
runs using two different numbers of particles per stream, 
Netr = 1.25 x 104 and Nety = 7.5 x 10°. We find that 
the inferred confidence intervals for all parameters are 
comparable in all cases. There’s perhaps a tendency for 
the confidence intervals obtained using a lower Ngt; to be 
10% wider, as a natural consequence of the compara- 
tively higher values of T;, which are a result of a slightly 
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Fic. 7.— Corner plots for the parameters recovered fitting the dog leg stream in NGC 1097: contours indicate the 68% and 95% confidence 


regions. 


higher Poisson noise in the effective pixels. 


5.1. The logarithmic total density slope 


Figure [7] shows the 68% and 95% confidence contours 
of the parameters that define the qualitative global mor- 
phology of the stream. As in the test runs, Fig. [6| the 
only noticeable degeneracy we identify is between the 
outer density slope y, and the break radius ro, through 
the orbital energy of the potential r..,/ro. The asymp- 
totic inner density slope is found to be 


4, =144£0.15, (15) 


where we are reporting median, 16% and 84% quantiles, 
identifying the 68% confidence interval. The stream’s 
morphology allows us to clearly see the density profile 
steepening from this central value, and the constraint on 
the asymptotic value of the outer slope is 


Yo =4.1+0.6. (16) 
The break radius is at 
SB Pes (17) 


well within the apocenter of the progenitor’s orbit and 
the total radial extension of the stream. 

The full radial profile we infer for the logarithmic to- 
tal density slope is displayed in Fig. |8| where both 68% 
and 95% confidence contours are shown. Vertical lines 
identify the inferred pericenter of the progenitor 


rea ioe (18) 
the break radius, and the orbital apocenter 


Pap = 1581, Roe; (19) 


The progenitor lies on a significantly eccentric orbit 
j = 0.12 + 0.08 , (20) 


and has just recently passed apocenter, 
t/(Tpp/2) = 5.23 + 0.07 , (21) 


where t = 0 here identifies the first pericentric passage 
and t/(T;»/2) = —1 is infall, or the first apocenter. It 
follows that, in the mean time, the progenitor has experi- 
enced three pericentric passages, at t/(T,»/2) € {0, 2,4} 

A sampling of the best-fitting model -A,¢ is displayed 
in Fig. [9] The left panel shows the position of the pro- 
genitor, in black, together with the distribution of stream 
members, within the same projected window probed 
by Fig. The stream’s morphology is easily recog- 
nised, with its characteristic X-shape and bright ‘dog 
leg’. Color coding indicates shedding time as from the 
displayed legend. For ease of reference, we have associ- 
ated a number to each of the four plumes of the stream. 
Plume 1 is the ‘leg’, contains the remnant of the pro- 
genitor, and is a combination of both leading a trailing 
material. It collects essentially all stars shed after the 
first pericentric passage. Plumes 2, 3 and 4 are all made 
of leading material lost at the first pericentric passage, 
and plume 4 is the farther away from the progenitor in 
terms of differential streaming, and — as a direct con- 
sequence — also the less extended. The right panel is a 
full map of the model residuals, normalised by the uncer- 
tainty of each effective pixel. As illustrated by the colour 
coding, residuals are everywhere quite low, rarely above 
3 in absolute value, and the only noticeable deficiencies 
of the model are: (i) a somewhat reduced brightness in 
plume 4 (visible in red in the residuals map), and (ii) a 
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Fic. 8.— The inference on the logarithmic density slope of the 
total density profile of NGC 1097. Grey shading identifies the 68% 
confidence region, shown together with the 95% confidence area. 
Vertical lines illustrate how the plotted radial interval compares 
to the inferred values of orbital pericenter, apocenter, and break 
radius of the density profile. For each, the 68% confidence region 
is displayed. 


mildly puffier 90 degree turn in plume 1 (evident as a 
blue clump). We return on this in Section 6. 

In order to illustrate in more detail how the stream’s 
projected morphology is capable of constraining the slope 
of the total density profile, we show the effect of varing 
the outer density slope y., while all other parameters 
are fixed as from the best-fitting model .4,¢. We explore 
approximately the edges of the 68% confidence region of 
‘Yo. As shown in the upper panes of Fig. [10] increasing the 
outer density slope to 7. = 4.7: (i) increases the speed 
of the differential streaming, so that plume 4 appears 
considerably brighter at the expenses of plume 2; (ii) 
implies significant increase in the galactocentric distances 
of the apocenters of the stream’s plumes, as for similar 
orbital energies larger distances can be reached thanks to 
amore slowly rising gravitational potential; (iii) increases 
the angles between plume 1 and each of the others, as a 
consequence of the change in the shape of orbits. The 
result is that the location and shape of all four plumes is 
changed enough to result in substantial residuals, shown 
in the upper residuals map with the same colour coding 
of Fig. [9] much higher than in the case of the best fitting 
model. 

The effect of a lower y, = 3.7 is analogous, although 
changes have the opposite qualitative direction, as proved 
by a comparison between the two residuals maps: the 
streaming speed is reduced and so are the length of the 
plumes and the angles between them. From this overview 
we see that the most important single pieces of informa- 
tion providing the constraints on the outer host density 
profile are: the length of plume 3 with respect to the 
current position of the remnant (or alternatively, to the 
length of plume 1), and the angles between plume 1 and 
all the others — it is worth recalling that the lengths of 
plumes 2 and 4 are not constrained by the data (see Sec- 
tion 2). 


5.2. The total mass profile 


The logarithmic density profile is a dimensionless quan- 
tity and can be derived from the surface brightness map 


of the stream alone. To associate a physical scale to the 
break radius of the profile, a dimensional scale of length is 
necessary and this is provided by the projected distance 
of the progenitor from the centre of NGC 1097, eqns. 
and (6). However, in order to proceed further, we need 
a dimensional scale to transform our inference on (7) in 
an inference on mass. In order to do so, we use the HI 
rotation curve measured by (2003): we 
make reference to the outer parts of the available rotation 
curve, which is observed to be approximately flat: 


V.(29.6 + 3 kpc) = 270+ 10 kms™' . (22) 


We use this condition to bring each of the dimension- 
less models accepted by our walkers to dimensional ones. 
Note that, in all inferences that require use of (22), re- 
ported uncertainties are the combination of the uncer- 
tainties deriving from the stream model, together with 
those propagated by the observational errors of the con- 
straint (22). which we assume have a Gaussian distribu- 
tion. 

Figure |11} shows the resulting circular velocity profile 
of NGC 1097, in both linear and logarithmic scale. The 
yellow point in the upper panel displays the normalisa- 
tion imposed by (23). As measured by[ Higdon & Wali 
(2003), the rotation curve is approximately Hat at those 
radii. The stream model naturally reproduces this find- 
ing and suggests the circular velocity remains flat out to 
r 2 ro. At even larger radii, however, the density profile 
steepens and the circular velocity declines considerably, 
as illustrated by the left panel. 

The combination of the constraint (22) and our stream 
model also provides a total mass: N 1097 is found to 
have a sizeable virial mass of 


logio(M200/Mo) = 12.23 + 0.12 , (23) 


and a virial radius of 
r200 = 2501?5 kpc . (24) 


We do not subtract stellar mass explicitly, and infer the 
concentration of this halo using the total mass profile. 
This implies 

C200 = Saar ‘ (25) 


Figure shows the degeneracy that exists between 
asymptotic outer density slope y. and concentration, and 
correspondingly, virial mass. Steeper density profiles re- 
quire larger break radii, implying lower values of the con- 
centration. In turn, less steep density profile can accom- 
modate larger values of the virial mass. 


6. RESULTS: THE PROGENITOR 
6.1. Mass evolution 


The main ingredients in the global morphology of the 
stream are the density profile of the host, which we have 
described in Section 5, and the mass of the progenitor. 
By themselves, the density profile and the progenitor’s 
orbit determine the approximate shape of the track of 
the stream. However, the actual length and width of 
the streaming material along this track are determined 
in first instance by the progenitor’s mass. Through the 
tidal radius (9), this fixes the speed of the differential 
streaming (see A15), and therefore the extension of the 
tidal tails with respect to the progenitor. 
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Fic. 9.— Left: a sampling of the best-fitting model, populated with Ngtr = 1.25 x 104 stream members and displayed within the projected 
region probed by the available data. The current position of the remnant of the progenitor is displayed in black. Right a map of the residuals 
of the best fitting model, colour coded according to the legend on the right. 


As mentioned in Section 3, we model the time depen- 
dence of the progenitor’s mass in an essentially non para- 
metric way, i.e. by using the masses at each pericentric 
passage as free parameters, and by interpolating with 
smooth splines between them. The progenitor is found to 
have recently passed the orbital apocenter (see eqn. (21), 
and has experienced three pericentric passages since in- 
fall. Using eqn. (22), we find that the apocenter preced- 
ing the first pericenter happened at 


ting = —5.4 = £0.6 Gyr , (26) 


and we can use this as an estimate of the infall time. 

Since that time, the total mass of the progenitor has 
evolved considerably: Figure shows the complete 
marginalised inference we obtain for the mass evolution 
of the progenitor. The upper panel displays the proba- 
bility distributions of the total masses at each pericen- 
tric passage. For construction, these are the masses con- 
tained within the tidal radius at that time, the probably 
distribution of which is shown in the intermediate panel. 
At its first interaction with the host, the progenitor is 
found to have a total mass of 


logo [Mperi1 (3-44 5'9 kpc)/Mo]| = 10.35 + 0.25. (27) 


Such a high mass is required in order to ensure the signif- 
icant differential streaming of the material shed around 
the first pericenter. As illustrated by the left panel of 
Fig. [9] such material entirely makes up plumes 2, 3 and 


At the second pericentric passage, we find that the to- 
tal mass is reduced by more than 1.5 orders of magnitude: 


logio [Mperi2(0.8619'33 kpc)/Mo] = 8.5193. (28) 


We interpret this as due to the eccentric progenitor’s or- 
bit: tidal forces likely strip very efficiently the biggest 
part of the progenitor’s dark matter halo when this 
plunges deep into the central regions of its host for the 
first time. Luminous material lost around the second 
pericenter is in the leg of the stream, as shown in Fig. 
and actually makes up both the peculiar ‘foot’ and putty 
‘knee’. After the first pericentric passage, the surviving 
remnant is compact enough to be quite resilient to tides, 
and mass evolution slows down considerably, with the to- 
tal mass at the third and last pericentric passage being 
just slightly lower than at the second (see Fig.|13). Ma- 
terial lost since the second pericenter has a much lower 
differential streaming and has not gained much distance 
from the progenitor, making up for the straight part of 
the plume 1, or the ‘leg’ itself. 

Simplifying, these constraints are provided by the dis- 
tances between progenitor and material shed at each 
pericentric passage, together with the extent of such ma- 
terial along the track of the stream. Both quantities are 
directly related to the speed of the differential streaming, 
which is a function of the tidal radius at the time of es- 
cape, and then a function of progenitor mass. It follows 
that uncertainties on the total mass of the remnant at the 
current time are comparatively larger, as the amount of 
stream members shed between the last apocenter and the 
current time is comparatively much smaller: 


logio [now (0.36456 kpc)/Mo] = 7.4798. (29) 


The complete constraints on the mass evolution of the 
progenitor are collected in the bottom panel of Fig. |13} 
which shows the total mass m at each of the three peri- 
centric passages, and today. Here time is measured from 
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the current time ¢ = 0, and the 68% confidence region 
on each pericentric passage is also shown. 


6.2. The ‘missing’ trailing tail 


As to the fractions of luminous material shed at each 
pericentric passage, it is unfortunately difficult to trans- 
form the inferences of the stream model into dimensional 
units. This would require resolved colour information 
along the stream. On a qualitative level, our modelling 
approach measures that the fractions of material shed 
at each pericenter decrease with time, and that the first 
interaction with the host strips the progenitor of a signif- 
icant amount of luminous mass (~ 0.6 of the luminance 
of the total shed mass). 

Note however, that while material shed at the second 
and third pericentric passages is entirely probed by the 
available data (see Fig. [9}, visible material shed around 
pericenter 1 is limited to the leading tail only. The cor- 
responding trailing tail lies entirely outside the region 
probed by Fig. Its location is illustrated in Fig. 
the trailing tail extends to over 100 kpc in projection 
and, as shown in the lower panel, its actual radial dis- 
tance from the centre of NGC 1097 reaches over 230 kpc, 
almost out to the nominal virial radius ro99. Our model 
assumes that equal amounts of material are lost from 
leading and trailing tidal radii. As shown by F is. [14 un- 
der this assumption the trailing tail contains a significant 
amount of mass and should be almost as bright as the 
other plumes. 

There are however at least two possible reasons as to 
why the trailing tail might have escaped detection so far. 
The first is that the hypothesis of equal escape fractions 
at the leading and trailing tidal radii is wrong for suffi- 
ciently close pericentric passages. Material could be lost 
in larger amounts from the leading tidal radius, where 
tides are stronger, so that our working assumption may 
well be over simplistic. The second possibility has to 
do with the exact profile of the density slope at large 
radii. Our model can account for significantly steep pro- 
files, but does not allow for truncations at even larger 
radii. A density slope that continues to steepen with 
radius would cause an additional increase in the radial 
distances reached by the trailing tail, up to determining 
the escape of part of such material. 


6.3. Line-of-sight kinematics 


Figure[14]displays a map of the line-of-sight velocity of 
the stream material, with respect to the systematic ve- 
locity of NGC 1097. Note that the models are symmetric 
with respect to the transformation 


Llos —* —DLos 5 (30) 


and additional information is required to fix the sign of 
the line-of-sight velocities, as well as the sign of the line- 


of-sight position of the progenitor at the current time 
eles (2010) mens that the progen- 
itor is receding faster than NGC 1097, and we fix the 
direction of the line-of-sight accordingly. In more detail, 


(2010) measure 


vers ~~ —30+ 30 kms’ , (31) 
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Fic. 10.— An illustration of the qualitative effects on the global 
morphology of the stream associated to variations in the outer den- 
sity slope yo. All other parameters are fixed according to the best 
fitting model (Fig.|9). Upper panels: the outer slope is steepened 
to Yo = 4.7, a corresponding stream and its residuals are shown. 
Lower panels: the outer slope is decreased to yo = 3.7. 


while the predictions of the model are 


Up,los = —514i? kms— 
Zpjlos = 1384417 kpe : 


1 


(32) 


so that the progenitor is nearer to us than NGC 1097. 
As shown in Fig. the entire leg is found receding, 
with a clear velocity gradient along its length. An even 
more marked gradient can be observed along plume 2, 
while plume 8 has a high internal velocity dispersion due 
to containing at the same time material that approaches 
and recedes our viewpoint. This is also clear from the 
bottom panel of the same Figure, which shows a plot 
of radial distance against line-of-sight velocity. Colour 
coding is the same as in Fig. |9} and shows the escape 
time of the stream particles. The extended trailing tail 
has mainly positive line-of-sight velocities, as plume 2; 
plume 3 instead is folded around its apocenter. 


6.4. Internal kinematics 


We find that the random velocities at escape increase 
with time, from o, + 1 kms~! for material lost at the first 
pericentric passage, to 0, © 3 kms~! at the current time. 
More interestingly, the model finds that non negligible 
ordered kick velocities are required to explain the mor- 
phology of the dog leg. These are larger than the random 
velocities and therefore correspond to internal rotation in 
the progenitor. Fig. [15] shows the ratio between the or- 
dered and random components of the kick velocities at 
escape for all stream members. The X-shape, defined by 
plumes 2, 8 and 4, is well described by material escaping 
with negligible coherent velocities. However, in order to 
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Fic. 11.— The circular velocity profile inferred for NGC 1097 in 
logarithmic scale (left) and linear scale (right). The yellow data- 
point in the right panel shows the observational constraint eqn 22}, 
used to normalise our dimensionless stream models. 


reproduce the dog leg in plume 1, values of OViot/os 2 1 
are required. 

Opposite signs of the coherent escape velocities for 
trailing and leading material allow them to stream to 
opposite sides of the main track of the leg. The trailing 
part of such material entirely makes up for the peculiar 
90 degree turn at the edge of plume 1, currently at apoc- 
enter (see also the bottom panel of Fig |14). Opposite 
kick velocities in the leading material cause the shift and 
peculiar morphology of the bright ‘knee’ of the dog leg. 
Both these features are the signature of internal rotation 
in the progenitor. The required rotation is prograde, and 
although limits on the precise inclination of the internal 
angular momentum are loose, we find that ding < —0.8, 
i.e. the disk is indeed somewhat inclined with respect to 
orbital plane. 

As noted in Section 5, our model seems to produce a 
dog leg that is slightly puffier than required by the data, 
and somewhat shorter, as shown by the residual map of 
Fig.[9| We advocate this is due to our certainly simplified 
treatment of the escape conditions, and in particular of 
the relation between escape velocities and spatial spreads 
in the initial conditions. The need for ordered kick veloc- 
ities increases the value of dv, which in turns makes the 
stream locally thicker because of the prescriptions (10). 

Finally, we can use our measurement of the remnant’s 
mass at the current time, eqn. to estimate the ve- 
locity dispersion of the remnant: 


gg lot ions * (33) 
where we have adopted the virial relation provided by 
(2011) and the half-light radius of the 
remnant as measured by (2010), r, + 300 


pe. 
7. DISCUSSION 
7.1. Limits of the model 


In this Section, we discuss some of the limitations of 
our modelling approach, and the possible biases they may 
introduce. 


7.1.1. Dynamical friction 


We start with dynamical friction, which is not taken 
into account by the model, so that the parameters of the 
progenitor’s orbit remain constant in time since infall. 
We first note that dynamical friction operates on the 
progenitor, but is negligible on the very light and diffuse 
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Fic. 12.— The degeneracy between the outer density slope yo and 


concentration c299, or virial mass M299. 68% and 95% confidence 
regions are displayed. 


material populating the stream once this is shed. For 
example, the morphology of the leading tail lost at the 
first pericenter — plumes 2, 3, and 4 — is unaffected by 
dynamical friction. As a consequence, the logarithmic 
slope of the total density profile of the host, which mainly 
derives from the shape of such tail, is largely unaffected 
by our working hypothesis. 

On the other hand, the orbital evolution of the pro- 
genitor driven by dynamical friction may cause shifts in 
the pericentric distances reached at each successive peri- 
centric passage. In turn, such shifts cause variations to 
the speed of the differential streaming of material shed 
around that time. In other words, while material lost at 
the second and third pericentric passages is assumed to 
be lost around the pericentric distance of the first pas- 
sage, its actual distance might have been smaller due to 
dynamical friction. This translates in a bias on the in- 
ferred masses of the remnant at the second and third 
passage. As the pericenter shrinks, the speed of the dif- 
ferential streaming increases (A15). As a consequence, 
masses that are even lower than in eqns. would 
be required after the first pericentric_passage — or con- 
versely an even higher initial mass 7. Therefore, dy- 
namical friction would make the mass evolution of the 
progenitor even starker. 

It turns out, however, that dynamical friction has lit- 
tle effect on the precise case in hand. We show this by 
providing an estimate of the relative loss of angular mo- 
mentum caused by each pericentric passage. We assume 


that the friction force can be described as in 
(1987), icc. as it were excerpted by a sea o 


Viaxwellian particles, with characteristic velocity v,. As 
the orbit of the progenitor is strongly radial, we approx- 
imate the interaction time with dt © 2rperi/vUp, where 
Up is the velocity of the progenitor at pericenter. Under 
these assumptions we find 


2 
oD oma bell (=) ; (34) 


Jp Us Tperi \ Up 


where A % rperivs/Gm is the usual Coulomb logarithm, 
and the precise value of the dimensionless constant c de- 
pends on the hypothesis made on the phase space distri- 
bution of the deflectors, c + 0.43 for a Maxwellian distri- 
bution. Even at the first pericentric passage, when the 
progenitor is at its most massive, we find a very limited 
effect, |OJp|/Jp S 0.05, comparable with the uncertainty 
of our measurement (see eqn. (20)). After the consider- 
able mass loss of the first interaction, dynamical friction 
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Fic. 13.— The mass evolution of the progenitor of the stream. 
Upper panel: marginalised pdfs for the total masses of the progen- 
itor within the tidal radius at the first (blue), second (green), and 
third (yellow) pericentric passages, together with the mass inferred 
for the remnant at the current time (red). Intermediate panel: with 
the same colour coding, the marginalised pdfs for the correspond- 
ing tidal radii. Lower panel: the total mass evolution in physical 
time. 


becomes negligible, with |dJp|/Jp S 0.002 at the second 
pericentric passage, and accordingly smaller thereafter. 


7.1.2. Host mass evolution 


Another limitation is connected to the assumption of 
a static host potential. If the mass M(r) increases as a 
function of time, this implies an evolution in the orbital 
energy of each streaming particle. For a stream like the 
one studied here, this effect is expected to be small. In a 
Milky Way sized halo, the mass within 100 kpc is already 


in place to within a few percent since redshift z = 1 
(Diemand et al.|/2007} |Boylan-Kolchin et al.|}2010). For 
our stream, infall is at z = 0.5 implying, on average, 


a very small mass evolution within the extension of the 
observed plumes. 

More in detail, a first effect of mass evolution is to 
make streams somewhat longer and thicker. The precise 
change in the energy of each particle depends on its orbit 
A, which couples it to the time dependence of the host 
mass evolution: 

OCF mney. aif OUD ay . (35) 
r,t) t 


Close stream particles with similar orbits experience a 
similar energy evolution, but for a stream that is long 
enough to have at the same time members in opposite 
orbital phases, eqn. introduces an additional energy 
spread. The energy spread of the shed material fixes 
the speed of the stream growth (A15), so that, for the 
same values of the progenitor mass and internal velocity 
dispersion, mass evolution can make long enough streams 
even longer and thicker. As a consequence, this could 
potentially cause an overestimate in the masses of the 
progenitor and in the random velocities at escape. 


The relevant dimensionless ratio that quantifies the 
importance of such effect is therefore the ratio between 
the additional energy spread introduced by mass evo- 
lution, o(F)m.ev., and the characteristic energy differ- 
ence (between tail and progenitor) that fixes the aver- 
age tidal streaming speed 6F;. For the particular case 
in hand, the orbital phases of stars lost after the first 
pericenter have not grown different enough from each 
other to result in an appreciable energy difference from 
mass evolution. Therefore, the only part of the stream 
that might be affected is the long leading tail lost at 
the first interaction. Even if we largely overestimate 
the spread in energy introduced by mass evolution with 
0(E) m.ev./E , © 6M/M, we obtain that the effect on the 
angular length of the stream is still very small, of the 
same order of 6M/M, i.e. a few percent. 

A second effect of mass evolution has to do with the 
shape of the orbits: if M(r) increases maintaining ap- 
proximately its spherical symmetry, the orbital angular 
momentum is conserved while energy decreases. This cir- 
cularises the orbits and determines a reduction of their 
apocentric distances with time, potentially affecting our 
measurement of the logarithmic density slope. This 
mainly derives from the ratios between the apocentric 
distances reached by the different plumes, so that when 
such reduction is uniform over all particles, no bias is 
created. 

Instead, cold dark matter halos evolve in an inside-out 
manner: at recent times their mass has evolved more 
slowly in their inner regions then in their outskirts (e.g. 
and references therein). As a result, it 
is the apocenters of the most energetic particles that are 
the most affected by mass evolution. While the potential 
well experienced by particles more deeply embedded in 
the potential stays unchanged, the particles with higher 
energies experience a reduction of their apocenters. The 
qualitative net effect is therefore a bias towards measur- 
ing lower values of y,: for a fixed energy gradient along 
the stream, a less steep density profile can compensate 
for the smaller differences between the apocenters of the 
different plumes determined by mass evolution. 


7.1.3. Spherical symmetry 


Finally, a more worrying limitation is that the present 
model uses spherically symmetric host density distribu- 
tions. Dark matter haloes have been shown to be sub- 
stantially non-spherical, and haloes with virial masses 
similar to what we infer for NGC 1097 typically have 
axes ratios distributed around a/c = 0.6 and b/c ¥ 0.8 
(Schneider et al.||2012| and references therein). On the 
other hand, the density distribution of NGC 1097 is cer- 
tainly non spherical in its central regions, which are likely 
dominated by its stellar disk and bar. 

The main effects of a non spherical potential are of 
two different kinds. First, as the equipotential surfaces 
are non spherical, apocentric distances are not constant 
in time, and vary with orbital phase. Second, as parti- 
cles experience a net torque, their instantaneous orbital 
plane is subject to a time dependent precession. If the 
non sphericity of the total density distribution is only due 
to the baryons, the quadrupole moment at the radii of 
the plumes is small, and so are these effects: the equipo- 
tentials at those distances are almost spherical (see for 


example|Amorisco & Bertin}2010), and the precession of 
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the orbital plane is not big enough to become apparent 
in a stream that extends for only three orbital oscilla- 
tions. However, if the density profile of the halo itself 
is non-spherical out to large radii, then these effects are 
substantial, and it may appear surprising that we are in 
fact capable of describing the stream almost down to the 
noise using a spherical model. 

We have made tests using single orbits, similar to 
those described in Section 4, in which mock datasets 
have been constructed by solving the equations of mo- 
tion for tracer particles in homeoidal density distribu- 
tions (Chandrasekhar 1960) These have initial condi- 
tions fixed by our inferences for the current position and 
velocity of the remnant, and integrate its orbit backward 
and foreword in time to mimic the four observed plumes. 
The trial density distributions we use are less prolate 
than suggested by cosmological simulations, with axes 
ratios a/c © 0.85 and b/c © 0.9, constant with radius, 
and oriented at random with respect to the present, in- 
stantaneous orbital plane. We find that use of our spheri- 
cal library to fit for such orbits is uniformly unsuccessful. 
Over three full orbital oscillations, orbits in these triaxial 
potentials loose their planar nature to the point that re- 
producing all four plumes (in this case apocenters) with a 
planar orbit from our library is impossible. In particular, 
we can never achieve y? levels that are comparable with 
those obtained in either our single orbit tests in spherical 
potentials, or in our stream fit. 

Although more detailed explorations are required, this 
suggests that, outside its central regions, the total den- 
sity distribution of NGC 1097 is not substantially non 
spherical. Observational constraints on the shape of dark 
matter haloes are quite loose, and almost completely 


limited to modelling the Sagittarius stream around the 
Milky Way’ (4 
and references therein). More recently, signif 
icant expectation is building up for the inferences ob- 
tained modelling the thin streams of globular clusters 
Bowden etal 2015), However, i 
comparative terms, these are limited to inner regions, 
where baryons are still important. On the other hand, 
it is widely accepted that dissipative dynamics and the 
lifecycle of baryons makes haloes rounder in cosmologi- 
cal hydrodynamical simulations (e.g. 
015) when 
compared to dark-matter-only runs. Our result provides 
observational evidence in this direction. 


7.2. Autopsy of a massive halo 


With a virial mass of ~ 2 x 10!?Mo, we find NGC 
1097 is not considerably more massive than the Milky 
Way. CDM expectations for the concentration of such a 
halo are c = 8'3°°, with the mentioned values indicating 


the 25% and 75% quantiles (Ludlow et al. [2014} 
STAN ULF NGC 1097 is in very good agreement wit 
these 


gures. In the stream’s morphology, we can see 
the total density profile bending from an inner slope of 
¥(Tperi) 1.5 around the pericenter, rperi * 4 kpc, and 
can clearly resolve its steepening. Interestingly, the cir- 
cular velocity remains approximately flat at radii that 
are even larger than the break radius of the total density 
profile, ro + 37 kpc. The very extended stream allows us 
to map the profile of the total logarithmic slope at least 
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Fic. 14.— Upper panel: a map of projected kinematics of the 
stream. Lower panel: line-of-sight velocity against galactocentric 
distance; as in Fig. ]9| colour coding indicates the shedding time, 
to facilitate the identification of the different plumes. 


until the orbital apocenter, rap, ~ 150 kpc, where we 
find y(rapo) = 3.9+ 0.5 (68% confidence region). To our 
knowledge, this represents the first statistical measure- 
ment of the outer density slope of a single extragalactic 
halo, i.e. on data that do not derive from stacking. 

As to the rather steep outer slope, our result certainly 
calls for the analysis of a wider sample of galaxies, so to 
establish whether NGC 1097 is just a case in the distribu- 
tion of outer slopes, or whether the outer density profile 
of galactic dark matter haloes is systematically steeper 
than expected. We recall that cosmological simulations 
show that the outer regions of dark matter haloes should 
be characterised by substantial halo-to-halo variations. 
Avila-Reese et al. suggest that the outer slope 
of galaxy-sized haloes correlates with environment and 
even when selecting isolated haloes only, they measure 
that the 68% confidence interval for the outer slope is 


sizeable, 2.5 < yo < 3.8 (see also|Prada et al.|2006). 
Recently, |Diemer & Kravtsov| (2014) have made clear 
that the simple NF'W profile is often a poor fit to the 


outskirts of ACDM haloes, even in the region 0.56 S 
r/T200.m <1. They show the significant halo-to-halo 
scatter in the outer density slope at z = 0, with many 
galaxy-sized haloes having envelopes that extend all the 
way down to ¥(1200,m) * —5. Their analysis concludes 
that the value of the density slope around the virial ra- 
dius is systematically associated to the mass accretion 
rate of the halo: all across the mass scale, steeper outer 
density profiles correspond to higher rates of recent phys- 
ical mass accretion. However, note that our apocen- 
tric distance is not large enough to explore the radial 
eng suidied in (2014), 150 kpc 
= 0.35 r200,m, underlying the need for a statistics of ex- 
tended streams and measurements of more outer density 
slopes. 
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Fic. 15.— The ratio between the ordered 6V,o¢, and random os, 
components of the escape velocity of each stream member (colour 
coding is illustrated by the legend). Both 90 degree turn and puffy 
‘knee’ in plume 1, the ‘dog leg’, are the signature of the internal 
angular momentum of the progenitor. 

In fact, it is quite surprising that using a simple, spher- 
ical broken power law profile we can describe the stream 
almost down to the noise. Apparently, the dynamics of 
the streaming material is not strongly affected by the 
details and local asperities of the gravitational potential, 
and a simple description like the one used here is suffi- 
cient to capture its salient properties. At the same time, 
our study shows that a full exploration of the parameter 
space is absolutely necessary: (2003) 
concluded that X-shaped streams could not be created 
within spherically symmetric hosts because their suite of 
N-body simulations did not produce any. The approach 
of this paper is considerably less numerically expensive 
than full N-body simulations, and therefore allows for a 
more systematic inspection of the wide range of param- 
eters. This is very promising and will allow the detailed 
study of more massive galaxies through the tidal features 
that encircle them. It certainly encourages the search for 
more low surface brightness, thin and extended streams. 
Here, we have made use of the projected morphology of 
the stream only, but kinematics of discrete tracers will 
be particularly helpful. 


7.3. The life of a cannibalised dwarf galaxy 


Galianni et al.| (2010) report that the remnant of the 
progenitor is highly symmetric, and has a well devel- 
oped core-halo structure. Its nucleus is bright (= 5.9 x 
10°Ly.o) and compact, but its half-light radius, certainly 
smaller than 80 pc, could not be measured precisely be- 
cause of seeing. The halo of the remnant, instead, is 
observed to extend for the full width of the stream, for 


2014). In the case of N 


a luminosity of ~ 2.6 x 10°Zy.o and an half-light ra- 
dius of a few hundred parsecs. Its spectrum suggests its 
metallicity is not too low, with [Fe/H] ~ —0.7 providing 
a reasonable fit. 

Our analysis reconstructs that this is the remnant 
of a disky dwarf, and that the progenitor was con- 
siderably more massive at infall, with a total mass of 
=~ 2.5 x 10'!°Mo within a few kpc (see eqn. (27). Infall 
happened ~ 5.4 Gyr ago, and in the mean time, the pro- 
genitor has experienced three pericentric passages. The 
first has been the most destructive, reducing its total 
mass by over 1.5 orders of magnitude. Most likely, this 
corresponds to the removal of the largest part of its ex- 
tended dark matter halo, leaving the remnant with only 
a few percent of its initial mass after just one interac- 
tion. The progenitor was internally rotating, as testified 
by the peculiar morphology of the dog leg, which requires 
ordered kick velocities 6V,o/o0, > 1 for the material lost 
around the second pericentric passage. 

Although estimating the total luminosity of the stream 
is not an easy task, from its approximate shape and mean 
surface brightness, obtain Lstr 
5x 10’Ly.o. This implies that, during its post-infall evo- 
lution, the progenitor has also lost the largest part of its 
stellar component. With a lower bound to the total stel- 
lar mass of © 10°Mo, and assuming that the virial mass 
of the progenitor at infall is at least Mj, © 10''Mo, such 
a dwarf galaxy is then roughly comparable with the ex- 
pectation of most abundance matching recipes (e.g. 


(et al 2010 Behroozi et al./2013). Furthermore, our pro- 


genitor is also in rough agreement with the relation be- 
tween nuclear stellar mass and total galaxy mass within 
the star’s characteristic radius. The relation constructed 
by using the galaxy population 
of the Virgo Cluster implies that a nuclear stellar mass 
of 2 x 10°Mo should correspond to a progenitor total 
mass of © 2 x 10° within the effective radius of its stellar 
distribution. 

Comparison between our measurement of the total 
mass of the remnant, eqn. (29), and the estimate of its 
stellar mass of © 6 x 10°Mo obtained by 
(2010) suggests that it may still contain some residua. 
dark matter in the central regions. The current remnant 
is not massive enough to be subject to substantial tidal 
friction (see Section 7.1), and will keep evolving almost 
undisturbed within the halo of NGC 1097. As ry, = rn, 
its diffuse halo will be entirely lost, most probably within 
a few Gyrs, and the only remains of the progenitor will 
be its bare compact nucleus. 

If we use the the current luminosity of the nucleus and 
the same mass-to-light ratio inferred by 
(2010) for the remnant, we get that the Roche lobe of the 
nucleus has a size of about a hundred parsecs when at 
the orbital pericenter. This implies that it will survive 
for a long time, details depending on the precise value 
of its poorly constrained half-light radius, r, < 80 pce. 
Tidal stripping of galaxy nuclei has been shown to be one 
of the channels of formation of ‘bright’ Ultra Compact 
Divan! ealatios 

UY, we are directly assisting 
to one of the stages in the formation of a low-luminosity 
UCD Norris et al 2014}, so that the 


progenitor of the deg leg stream proves that such evo- 
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lutionary path also determines the formation of a low 
luminosity tail in the population of non-globular stellar 
clusters. This study shows that analysis of the stream 
properties is capable of providing a one-to-one associ- 
ation between the observed remnant and the structure 
and nature of the progenitor dwarf galaxy, through the 
reconstruction of its post-infall evolution. 


8. SUMMARY 


By modelling the two-dimensional surface brightness 
map of its stellar stream, we have reconstructed the ac- 
cretion of a dwarf galaxy on to the massive disk galaxy 
NGC 1097. We have shown that modelling the projected 
map of the stream alone allows us to constrain at the 
same time the initial properties of the progenitor galaxy, 
its post-infall evolution, and the density profile of the 
host. We find that: 


e We can fully describe the properties of the stream 
using a spherical model for the total density distri- 
bution of the host galaxy. We perform tests using 
single orbits in triaxial potentials and find that this 
would not be possible if the halo of NGC 1097 were 
to be substantially triaxial out to large radii (i.e. 
with constant axis-ratios). 


e The progenitor dwarf galaxy is on a strongly ra- 
dial orbit, so that the stream extends out to over 
150 kpc, allowing us to map the logarithmic slope 
of the total density profile y of the host out to 
= 0.6r299. The shape of the stream implies that 
the density slope bends from an inner value of 
Y(Tperi * 4 kpc) = 1.5+0.15, to a quite steep 
value of 7(Tapo © 150 kpc) = 3.9 + 0.5. 


e We can clearly resolve this steepening with ra- 
dius and find that the location of the break ra- 
dius is in very good agreement with the predic- 
tions of ACDM, so that NGC 1097 has a concen- 
tration of cy99 = 6.773 and a virial (total) mass 


of Maoo = 1.8703 x 10!27Mo. 


e Our model predicts that the remnant of the progen- 
itor has a line-of-sight velocity of vpjos = -515i4 
kms! with respect to NGC 1097, which is in ex- 
tremely good agreement with the measured value 
of oe 


p,los — —30 + 30 kms7! (Galianni et al. 2010). 


e The progenitor dwarf galaxy, which is now observed 
as a nucleated dwarf with a total stellar luminosity 
of © 3.3 x 10’Ly,o (Galianni et al.|/2010), has in 
fact lost more than two orders of magnitude in total 
mass since infall, 5.4+0.6 Gyr ago, and experienced 
three pericentric passages. Our inference for its 
initial mass is log )|[m(< 3.4+1kpc)/Mo] = 10.354 
0.25. 


e The peculiar morphology of the dog leg stream 
is due to the initial properties of such progenitor 
dwarf, which had a non negligible internal angular 
momentum and was then likely a disky dwarf. Pro- 
grade rotation is mirrored in the coherent escape 
velocities that are necessary for reproducing both 
‘foot’ (i.e. sharp 90 degree turn) and ‘puffy knee’ 
in the dog-leg morphology, material lost around the 
second pericentric passage. 


e The remnant today is observed to have a marked 
‘core-halo’ structure. Our inference for the current 
total mass in the remnant implies that the stellar 
‘halo’ is not compact enough to be resilient to tides 
and will be lost entirely. Conversely, the bright 
and compact nucleus is unaffected by both tides 
and dynamical friction and will survive for a long 
time. 


e With a luminosity of 5.9 x 10°Ly.o, our model 
shows explicitly that tidal stripping of dwarf galax- 
ies down to their nuclei is one of the mechanisms 
that brings to the formation of UCDs, in this case 
of what would be called a low-luminosity UCD, and 
can reconstruct the properties of its progenitor. 


These results encourage the observational effort to- 
wards the search for more extended, reasonably thin, 
stellar streams around massive galaxies. The combina- 
tion of a wider sample of systems and the statistical mod- 
elling approach presented here will allow us to unveil the 
properties of both massive haloes and dwarf galaxies. 


It is a pleasure to thank Aaron Romanowsky for con- 
structive comments. The DARK Cosmology Centre is 
funded by the DNRF. 
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